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The extraction of the strangeness form factors from parity violating elastic electron-proton scat- 
tering is sensitive to the electromagnetic form factors at low Q 2 . We provide parameterizations 
for the form factors and uncertainties, including the effects of two-photon exchange corrections to 
the extracted EM form factors. We study effect of the correlations between different form factors, 
in particular as they impact the parity violating asymmetry and the extraction of the strangeness 
form factors. We provide a prescription to extract the strangeness form factors from the asymmetry 
that provides an excellent approximation of the full two-photon correction. The corrected form 
factors are also appropriate as input for other low-Q analyses, although the effects of correlations 
and two-photon exchange corrections may be different. 

PACS numbers: 



I. INTRODUCTION 

The parity-violating (PV) asymmetry in elastic scat- 
tering of polarized electrons from unpolarized protons 
can be used to extract information on the strangeness 
contribution to the proton form factors [H, 0, HI- Be- 
cause the electromagnetic (EM) coupling is proportional 
to the quark charge-squared, scattering from the proton 
is strongly dominated by interaction with the up quarks. 
Electron-neutron scattering provides a different relative 
weighting of the up and down quark distributions, al- 
lowing one to study the difference between up and down 
quark contributions to the nucleon form factor under the 
assumption that the up-quark distribution in the proton 
is identical to the down quark distribution in the neu- 
tron, and neglecting heavier quarks. Because the parity 
violating cross section comes from interference between 
photon and Z exchange, the quark flavors have a differ- 
ent weighting in the interaction, allowing separation of 
up, down, and strange contributions to the form factors 
by combining proton and neutron electromagnetic form 
factors and parity-violating e-p scattering. However, the 
small contribution of the strange quarks to the parity- 
violating asymmetry requires precise knowledge of the 
contributions from the up and down quarks before one is 
able to achieve sensitivity to the strange quark contribu- 
tions. 

The parity-violating asymmetry arises due to interfer- 
ence between photon exchange and Z exchange, and in 
the Born approximation is given by Q 
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where G F is the Fermi constant, a is the fine structure 
constant, and Q 2 is the four-momentum transfer squared. 
The individual asymmetry terms can be written in terms 



of the proton's EM vector form factors, Ge p and Gmp, 
and the proton's neutral weak vector and axial form fac- 
tors, 
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where 9 e is the electron scattering angle, r = Q 2 / (4M 2 ), 
t _1 = (l + 2(l + r)tan 2 9 e ), 6w is the weak mixing angle, 
and e' = y/r(l +t)(1 -e 2 ). 

In the standard model and with the assumption of 
isospin symmetry, the weak form factors can be expressed 
in terms of the proton and neutron EM form factors and 
the strangeness contribution to the nucleon EM form fac- 
tors, Ges and Gms, neglecting contributions from heav- 
ier quarks Making this substitution, and removing 
the common factor A = — {GeQ 2 )/{^o.^/2), Apy con- 
tains the terms 
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which depend only on quantities that are measured or 
which can be reliably estimated, and one final term, 
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which contains the unknown quantities of interest: Ge s 
and Gms- In the above expressions, we have written the 
denominator in terms of the e-p reduced cross section, 
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To extract the strangeness-containing term, the best 
known values and uncertainties for the other terms are 
needed. We present an analysis of the world e-p and e- 
n scattering data to determine the nucleon form factors, 
the e-p reduced cross section (the denominator of Eqs.[3} 
[6]), and their uncertainties. We also study the impact of 
the correlations between the different form factors as well 
as the effect of two-photon exchange (TPE) corrections 
on the extraction of Ge s and Gms- Electroweak radia- 
tive corrections have been calculated 0, @] , and their un- 
certainties do not generally limit the extraction of the 
strangeness contributions. 

The extracted form factors and uncertainties are also 
appropriate for use in the analysis of other high precision, 
low Q 2 experiments. However, the analysis of correla- 
tions and TPE exchange effects presented in this paper 
is aimed specifically at parity-violating clastic electron- 
proton scattering. Care must be taken in using these fits 
in analysis of other experiments, as it is necessary to de- 
termine if the analysis requires the Born form factors or 
simply needs the form factors as a parameterization of the 
elastic cross section. When the cross section is required, 
using the Born form factors requires making an explicit 
correction for TPE effects 0). For other cases, such as 
the extraction of the axial form factor from neutrino scat- 
tering Q, determining corrections to hyperfine splitting 
in hydrogen Q , or determining the Bethe-Heitler term in 
the analysis of DVCS measurements, one needs to care- 
fully consider whether TPE corrections are needed and 
to which degree they are different from the ones needed 
for the unpolarizcd cross section. 



II. ANALYSIS OF LOW-Q DATA 

A. Proton form factors 

A fit to the world e-p cross section data at very low 
momentum transfer has been described in [l(|. This fit 
uses a Continued Fraction (CF) expansion, 
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of Ge p and Gmp most suitable for the lower momentum 
transfers, and extends up to Q = y/Q 2 ~ 0.8 GeV/c. 
Note that these fits should only be used in the quoted 
range of Q values. The analysis includes the effects of 
Coulomb distortion which, contrary to common belief, 
are not negligible [ll[ . The effect of two-photon exchange 
beyond Coulomb distortion, which includes only the ex- 
change of an additional soft photon, has also been stud- 
ied fl2l . [l3j . In our main analysis, we will correct the 
proton cross sections for Coulomb distortion, though we 
also provide parameterizations using the full calculation 
for TPE effects and discuss the impact of TPE correc- 
tions on the neutron form factors. 



TABLE I: Fit parameters for the low-Q form factors, valid up 
to Q — 1 GeV/c, using the CF parametrization of Eq. [7] (with 
Q 2 in (GeV/c) 2 ) for Ge p , Gm p /h p , GMn/Hn, and the param- 
eterization of Eq. [5] for Geu- The proton data are corrected 
for Coulomb distortion. 





bi 


b 2 


63 


b 4 


h 


Gep 


3.440 


-0.178 


-1.212 


1.176 


-0.284 


G]\ip/ 


3.173 


-0.314 


-1.165 


5.619 


-1.087 


Geu 


0.977 


-20.82 


22.02 






Gmtl/ )J"n 


3.297 


-0.258 


0.001 







Here, we extend this fit to higher momentum transfers, 
up to Q = 1.2 GeV/c, such as to sufficiently bracket the 
kinematics covered by the different PV experiments. The 
approach taken is identical to the one employed in [lol ]: 
We start from the world cross sections for e-p scatter- 
ing El JSJS E3, EE El, M, M, El m, M, S 

(27|, [23, 29, [301, apply the Coulomb corrections accord- 
ing to ll|, |3l| and fit the cross sections with CF pa- 
rameterizations of both Gep and Gmp- The longitudi- 
nal/transverse (L/T)-separation is done implicitly by the 
fit. 

This approach allows one to keep track of the random 
and systematic uncertainties of the data and propagate 
them to the final quantities of interest. From the result- 
ing CF parameters and the error matrix of the fit one 
can calculate the values of Ge p and Gmp together with 
their random error for any desired value of Q. To obtain 
the systematic error of the proton form factors, each data 
set is changed by the quoted systematic uncertainty, the 
world set is refit, and the change in Ge p and Gmp cal- 
culated. These changes are added up quadratically for 
all data sets, yielding the systematic uncertainty on Ge p 
and Gmp- This is usually the dominant error. The to- 
tal error is obtained by adding quadratically the random 
and the systematic errors, a procedure that should be 
applicable given the large number of data sets used. 
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FIG. 1: Uncertainties in the fits for Ge p (solid) and Gm p /Mp 
(dashed). The uncertainty is the random and systematic un- 
certainties, combined in quadrature. 
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The fit to the proton cross sections, after correcting 
for Coulomb distortion, yields the coefficients given in 
Table H The fits are valid for Q from 0.3 to 1.0 GcV/c. 
The uncertainty for Ge p {G m p / Vp) , is given by the solid 
(dashed) line in Fig.[TJ At low Q, the error bar on Gmp 
is larger than the one on Gep as the data base is less 
complete, although in the low-Q region there are two 
data sets with measurements at 180° [22l. |26I|. 



B. Neutron form factors 

The most precise data for the neutron magnetic form 
factor Gmu come from measurements of the ratio of 



2 H(e, 



en. 



to 2 H(e,e'p) [3J, |33j. The value of G Mn is 



extracted from the neutron cross section, which is de- 
termined from the combination of the neutron to proton 
ratio in deuterium and the (free) proton elastic cross sec- 
tion. We also include measurements from the asymmetry 
on polarized 3 He (34j . which are of somewhat lower pre- 
cision, and data points for Q < 1.3 GeV/c from Ref. [351 ] . 
The high Q points have larger uncertainties, and are out- 
side the range of validity of the fit, but are included to 
avoid "extreme" behavior for Q < 1 GeV/c. We fit these 
data to a 3rd order CF expansion, and the parameters 
are shown in Tab. [U The random and systematic uncer- 
tainties of Gmu have been estimated in Ref. [Hj]. In the 
range of momentum transfer of interest here, they are 
approximately 1.5%, roughly independent of Q. 

The value of the neutron charge form factor, Geh, 
is obtained by fitting all data presently available from 
polarization-transfer experiments [H, [33, [H, [H, 0, EH 
E2, EH, EH- Care has been taken to employ the most 
recent values, as some of the experimental Geu's pub- 
lished early on did not contain the best corrections for 
FSI and MEC (or no corrections at all). To study the un- 
certainty due to FSI and MEC corrections, an additional 
uncertainty equal to 30% of the calculated correction was 
added in quadrature with the experimental uncertainties. 
Including this additional uncertainty in the extraction of 
Geu has little effect on the fit or the uncertainties. Also 
included in the fit are the G En -v alues determined from 
the deuteron C2 form factor (45[ and the slope of Gb„ 
at Q = 0, known from n-e scattering. For Geu the er- 
ror bars of the published data contain a mix of random 
and systematic uncertainties. This mix is difficult to take 
apart, and therefore no distinction between random and 
systematic errors is made here. 

The fits of Geu are done using a modified 3-parameter 
CF expansion: 



GEn{Q) = 0.484 • Q 2 ■ Gcf, 



(8) 



with Q 2 in (GeV/c) 2 . The constant value in front fixes 
the slope at Q 2 = to match the measured rms-radius 
squared value of -0.113 fm 2 (46[. The fit parameters are 
given in Tab. (H and the fit is valid up to Q = 1 GeV/c. 

The error matrix is used to compute the error of Geu 
at any desired value of Q. Fits using a functional form 
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FIG. 2: Total uncertainty for the fit to Geu- 



similar to the Galster fit [47| are systematically below our 
fit, while fits attempting to include an explicit pion cloud 
contribution [48| lie above our fit. We performed fits 
using different fit functions and took the fit dependence 
as an additional contribution to the uncertainty in Geu- 
This is the dominant source of the uncertainty below Q ss 
0.3 GeV/c, where no direct measurements exist. The 
final estimated uncertainty on Geu is shown in Fig. [21 
The figure shows the absolute uncertainty on Geu] the 
relative uncertainty is well parameterized by taking the 
minimum of (5.2 + 12.6 -Q 2 )%, which fits the curve below 
Q = 0.3 GcV/c, and (9.2 + 11 ■ exp (-Q 2 /0.19))%, which 
fits above Q = 0.3 GeV/c. 



C. Cross Section 

In the Born approximation, Apy is given by Eq. [1] 
The inclusion of two-photon exchange terms leads to 
the replacement of the Born form factors Ge p (Q 2 ) and 
Gm p {Q 2 ) with generalized form factors that depend on 
both £ and Q 2 , as well as introducing two new terms, A' M 
and A' A Given a complete calculation of the two- 
photon exchange correction, one can extract the Born 
form factors by correcting the Roscnbluth and polariza- 
tion extractions for TPE effects, and then applying the 
TPE corrections to Apy- However, the TPE corrections 
to the denominator of Eq. Q] are identical to the correc- 
tions to the e-p unpolarized cross section. So rather than 
correcting the unpolarized cross section measurements 
for TPE and then re-applying TPE correction to eval- 
uate (J re d , one can make a model-independent evaluation 
of er rc d by taking a fit to the TPE-uncorrected e-p cross 
section. Thus, we also provide a fit to the measured e-p 
cross section, without applying any kind of TPE correc- 
tion. 

The procedure is identical to the extraction of the pro- 
ton form factors, except that a global fit is performed to 
the uncorrected cross sections. The reduced cross section 
is fit to the the form a ICI ± = TF m (Q 2 )+eF e (Q 2 ), such that 
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TABLE II: Fit parameters for the low-Q e-p cross section, 
neglecting TPE corrections. 
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b 2 


63 


64 


65 




3.366 


-0.189 


-1.263 


1.351 


-0.301 


F m 1 fip 


3.205 


-0.318 


-1.228 


5.619 


-1.116 



in the Born approximation, F m = Gm p and F e = Ge p - 
While the fit to the TPE-uncorrected data can also have 
an s dependence, a global analysis of the e dependence 
of a ep indicates that that deviations from linearity are 
extremely small [4!| . Table [TT| gives the parameters for 
the fit to the uncorrected cross sections. This fit is ap- 
propriate both for the reduced cross section term in the 
evaluation of Apy, but also as a parameterization of the 
clastic cross section with the TPE corrections absorbed 
into the fit function. It is therefore useful as a low Q 
model of the elastic cross section if one does not wish to 
explicitly treat the TPE corrections for unpolarized e—p 
scattering. 
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FIG. 3: Total uncertainty for the fit to the TPE-uncorrected 
value of cr re d at the kinematics of several past and planned 
measurements. The kinematics of forward angle are shown 
for the JLab (O), while backward angle kinematics are shown 
for Bates (□) measurements. The Mainz (+) points are for 
forward angle measurements at low Q, and backward angle 
measurements at higher Q. 

The uncertainties in F e and F m are essentially identical 
to those of Ge p and Gm p - However, the calculation of the 
uncertainty in cr ro d requires special care. The values of F e 
and F m are strongly correlated as they result from the 
(implicitly made) L/T-separation of the cross sections. 
The uncertainty on the cross section thus is smaller than 
the one one would obtain by combining the errors in F e 
and F m . The error matrix is used to evaluate the random 
uncertainties, while the systematic uncertainty is taken 
as the combined effect of individually varying the normal- 
ization of each data set. It is not possible to provide a 



TABLE III: Fit parameters for the low-Q proton form factors, 
using the two-photon exchange correction from Ref. P3] ■ 
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63 


64 
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Gep 
Gmp 1 


3.478 
3.224 


-0.140 
-0.313 


-1.311 
-0.868 


1.128 
4.278 


-0.233 
-1.102 



simple parameterization for the cross section uncertainty 
at all e and Q values. The cross section uncertainties at 
kinematics corresponding to a variety of PV experiments 
are given in Figure [3l For the forward angle measure- 
ments, the cross section uncertainty is a relatively simple 
function of Q, so the uncertainty for other large-e can be 
estimated from the small angle data in Fig. [3] The Q 
dependence is more complicated for large angle, and we 
have tried to include the complete set of planned mea- 
surements. 



D. Two-photon exchange beyond Coulomb 
distortion 

In the extractions of the proton form factors de- 
scribed above, we have applied Coulomb distortion cor- 
rections [ll|, EH to the e-p scattering data, and no cor- 
rection to the e—n data. Coulomb distortion takes into 
account the effect of a second soft photon, but does not 
include the contribution from a second hard photon. At 
low Q, the difference between Coulomb distortion and 
full TPE is small \Ej^. 

To gauge the effect of the full two-photon corrections, 
we perform another extraction of the proton form factors, 
after correcting a ep using the TPE calculation of Blundcn 
et al. fl2| , rather than the Coulomb distortion correction. 
These authors calculate the contribution of the exchange 
of a second photon, soft or hard, restricted to the case for 
the intermediary state being a proton in its ground state. 
This calculation includes only the unexcited proton in 
the intermediate state; the contribution from excited in- 
termediate states is neglected. This calculation explains 
most but not all of the discrepancy between Rosenbluth 
and Polarization measurements above Q 2 =2 (GeV/c) 2 , 
but appears to be sufficient at lower Q 2 values. A cal- 
culation including an intermediate A [51| also indicates 
that this contribution is important at Q 2 =2-3 (GeV/c) 2 , 
but provides only a small modification below 1 (GeV/c) 2 . 

Table |ni] shows the results of the fit to the proton cross 
sections corrected for TPE, as opposed to the Coulomb 
distortion corrections used for the fits in Table Q] For 
this fit, the explicit Coulomb distortion corrections done 
above were omitted, as the TPE corrections already con- 
tain the contribution from Coulomb distortion. The fit is 
valid for Q from 0.3 to 1.0 GeV/c. The uncertainties of 
Ge p and Gm p (random plus systematic added quadrat- 
ically) are essentially identical to the fit with Coulomb 
distortion (Fig. [T]) . While some of the parameters in the 
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fits are noticeably different, the difference between cor- 
recting for Coulomb distortion and TPE on the extracted 
form factors is small. While the size of the Coulomb dis- 
tortion corrections can be up to 3% for Ge p and 1% for 
Gmp, the difference between Coulomb and full TPE cor- 
rections is typically 0.3-0.4%, and never more than 0.5% 
for Ge p and 0.7% for Gm p . This difference is always 
less then the uncertainties shown in Fig. [TJ and more 
importantly, is well within the radiative correction un- 
certainties assumed in the initial measurements, which 
are dominated by the uncertainty in TPE contributions. 
While the estimates of the uncertainties in the radiative 
correction procedure were clearly underestimates when 
neglecting TPE corrections, they provide a reasonable 
estimate of the uncertainty in the TPE calculation, es- 
pecially for these small Q values. 

For the neutron, the TPE correction to the cross sec- 
tion as calculated by Blunden et al. fl2| is well param- 
eterized at low Q as Acr/er = 0.8% • Q 2 ■ (1 — e), with 
Q 2 in (GeV/c) 2 . For Q < 1 GeV/c,this yields a maxi- 
mum correction go Gmu of 0.4% at Q = 1 GeV/c and 
large scattering angle. Since most measurements are at 
relatively forward angle, the typical correction is sJO.1%. 

The calculated two-photon corrections to the polar- 
ization measurements of the neutron electric form factor 
are extremely small. It must be noted, however, that 
most modern experiments determining Geu measure an 
asymmetry depending on the ratio Geu/Gmu, so a TPE 
correction to Gmu propagates into the extracted value of 
Geu- As this correction is well below 1%, the effect is 
negligible compared to the experimental uncertainties of 
these measurements. 

Thus, the overall difference between the full TPE cor- 
rection and the Coulomb distortion is quite small. For 
protons, this different can amount to about half of the 
final uncertainty coming from the input form factors, but 
is within the radiative correction uncertainty applied in 
the individual measurements, and thus is properly ac- 
counted for in the final uncertainties. However, one also 
has to consider correlations between the uncertainties 
in different form factors, which can enhance the effect 
on the parity violating asymmetry. This is discussed in 
Sec. 



III. DETERMINATION OF A PV 

Given the nucleon electromagnetic form factors, as well 
as 9w and G\, one can calculate Ap = v ° 1 the PV asymme- 
try for Ges = Gm s = 0. One then takes the difference 
between the measured asymmetry and Apy° as the con- 
tribution from the unknown term (Eq. [5]). To obtain a 
reliable value for the strangeness contribution, one must 
include radiative corrections in evaluating Ap v ° , and de- 
termine the uncertainty in Ap^P due to uncertainty in the 
form factor and other terms. The uncertainty in Ap^P 
is usually determined by varying the individual form fac- 
tors that contribute to the asymmetry by their assumed 



uncertainties. This approach ignores two effects which 
could be significant in these measurements. First, it ne- 
glects the correlation between the extracted values of the 
electromagnetic form factors, which can impact the to- 
tal uncertainty on Apy°. Second, it neglects the impact 
of two-photon exchange corrections on Ap v °, as well as 
their effect on the extracted values of the form factors. 
In the following sections, we will study the effect of the 
correlated uncertainties between the extractions of the 
different form factors, and estimate the size of TPE cor- 
rections. We will present a procedure for determining 
the size and uncertainty of the parity-violating asymme- 
try that does not require an explicit calculation of TPE, 
but which minimizes the uncertainty in extracting the 
strangeness contributions. We will compare this to the 
result obtained if one ignores both TPE corrections to 
Apv and the TPE corrections in the extraction of the 
electromagnetic form factors, as has been done in all pre- 
vious extractions of the strangeness contributions. 

A. Impact of correlations 

The effects of these correlations need to be evaluated to 
obtain an accurate measure of the uncertainty in Ap v °. 
Taking the correlations into account can noticeably in- 
crease or decrease the contribution of the form factors 
to the total uncertainty. We examine here the impact of 
these correlations on the evaluation of Apy°, in the re- 
gion of 0.3 < Q < 1.0 GeV/c, where such measurements 
have been carried out or proposed. 

At the Q values of interest here the main issues of con- 
cern are the anti-correlation between Ge p and Gm p ex- 
tracted from Rosenbluth separation, the correlation be- 
tween Gmu and a ep (cr re d) for data extracted in measure- 
ments of the proton/neutron ratio or the 3 Hc quasiclastic 
asymmetry and the correlation between Geu and Gmu 
in polarization measurements that extract Geu/Gmu- In 
some cases, it is difficult to precisely quantify the level 
of correlation and difficult to propagate to the value of 
Apy°. The aim of this analysis is to determine where 
these correlations can be neglected or treated in some 
approximate fashion, and to determine the uncertainty 
related to these approximations. 



1. Correlation between Gep and Gm p 

At the low values of Q of interest here, the proton form 
factors are determined by L-T separations, which yield a 
significant anti-correlation between the extracted values 
of Ge p and Gm p . This has an impact in determining the 
uncertainty in er r cd; which appears in the denominator of 
most of the terms, as well as introducing a correlation 
between the errors in the term involving Ge p and the 
term involving Gm p - 

The largest effect results from the fact that the un- 
certainty on <7 re d is much smaller than one obtains by 
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varying Ge p and Gm p individually. By treating <7 re d 
and <5er ro d as being independent quantities from Ge p 
and Gmp, we eliminate the overestimate of the uncer- 
tainty, and we better account for TPE effects as well (see 
Sec. Ill Cp . When er re d is extracted directly from the cross 
sections, the remaining correlation between cr re d and the 
individual form factors is very small. 

The remaining effect is the correlation between the 
terms in Apy involving Ge p and Gm p (Eqs. [3] and [J). 
These terms have the opposite sign in the final asym- 
metry, so the anti-correlation between the values of Ge p 
and Gmp will tend to increase the total uncertainty. For 
small £ values, the term involving Ge p is only a few per- 
cent of the total asymmetry, and so its uncertainty has a 
negligible small effect (< 0.1% of Ap\f for e < 0.05). At 
large e values, this term contributes roughly 20% the of 
Apy , and so has a greater impact. If taken to be 100% 
anti-correlatcd with Gm p , the difference at large e grows 
from 0.1% of A s py° at Q 2 = 0.1 (GeV/c) 2 to 0.4% at 
1 (GeV/c) 2 . In reality, the effect will be smaller, as the 
correlation is not 100%, especially in the region where 
Gm p is mainly given by the data taken at 180°. The 
completed and proposed measurements of Ap^y° typically 
have ~10% precision on Apy , and never better than 4%, 
and so neglecting this correlation will again have a very 
small effect on the final uncertainty. Thus, it is a good 
approximation to treat Ge p and Gm p as uncorrelated, as 
long as one takes the uncertainty <5<7 ro d directly, rather 
than calculating <5er rc d from SGe p and SGm p - 



2. Correlation between Gmu and <7 r ed 

The most precise values of Gmu come from measure- 
ments of the ratio of d(e,e'n) to d(e,e'p). By normaliz- 
ing to a ep , measured on the proton, one can extract o~ en 
and thus Gmui since the contribution from Geu is almost 
negligible. Measurements utilizing quasielastic scattering 
of polarized electrons from polarized 3 Hc are essentially 
measuring the same quantity [34[ ■ The transverse asym- 
metry for scattering from the polarized neutron is nearly 
independent of the neutron form factors, and so the 3 He 
asymmetry is mainly sensitive to the dilution due to the 
two (nearly) unpolarized protons, and thus is sensitive 
to a ep /a en . Therefore, both experiments yield a direct 
correlation between the extracted value of Gmu and the 
value of o~ ep used in the analysis. 

Because of this correlation, it is important that TPE 
are treated in a consistent fashion. If the proton form 
factors are corrected for TPE, then the TPE contribu- 
tions must be included in calculating a ep as observed 
in the o- en /a ep measurements. Because TPE corrections 
were neglected in both the extraction of the proton form 
factors and the calculation of a ep as used in the Gmu ex- 
tractions, one obtains a correct parameterization of the 
unpolarized e-p cross section, as in Sec. Ill CI 

The typical uncertainties in <r ep at kinematics where 
Gmu has been extracted arc «1.4%, yielding a con- 



tribution to the uncertainty in Gmu °f sa0.7%. For 
Q < 1 GeV/c, the uncertainty in Apy due to Gun is 
close to half the size of the uncertainty coming from cr ro d, 
so a perfect correlation would have the effect of reducing 
this contribution to the uncertainty by roughly a factor of 
two. However, the uncertainty on er rc d is usually not the 
dominant contribution, and so the effect of reducing this 
contribution is never more than 0.4% of Apy . In fact, 
the effect is even smaller since the uncertainties are not 
100% correlated between different extractions of Gmu- 



3. Correlation between Ge-ti and Gmn 

The polarization measurements are sensitive only to 
the ratio GeuIGmu-, and thus the error in Gmu used in 
the analysis yields an identical shift in Geu- However, 
the contribution to SGeu from SG]\i n is a very small part 
(typically 10%) of the total uncertainty in Geu- 



B. Two-photon exchange corrections to Apy 

A full calculation of the TPE corrections to Apy re- 
quires starting with the TPE-corrected (i.e. Born) form 
factors, and then applying the full TPE corrections for 
parity-violating scattering to Eq. [T] We have made fits 
to the proton and neutron form factors, both with par- 
tial TPE corrections, neglecting the effect of a second 
hard photon, and with full TPE corrections. The full 
corrections are more model dependent, but are close to 
the corrections applied, as discussed in Sec. Ill Dl The 
uncertainties assumed for the radiative corrections are 
now consistent with the corrections applied, even for the 
partial correction. 

We apply TPE corrections based on the formalism by 
Afanasev and Carlson Q . This includes the effect of the 
two photon box (and crossed-box) diagrams, but not the 
effect of the 7-Z box diagram, which has been examined 
(for Q 2 = 0) in Ref. [52|]. It should be noted that for 
these corrections, it is not sufficient to apply the correc- 
tion for only the second soft photon; one must go beyond 
Coulomb distortion. This is because the Coulomb distor- 
tion is a long range contribution that, to first order, yields 
a helicity-independent rescaling of the cross sections, and 
thus cancels in the evaluation of Apy. So for Apy, one 
must use the full calculation, including the exchange of 
a hard photon. 

For convenience, we separate the TPE effects into three 
categories. First, there are two new terms, A' M and A' A , 
that appear in the expression for Apy. Second, the Born 
form factors that go into the terms Ae, Am, and A a 
in Eq. Q] are replaced with generalized form factors that 
depend on both e and Q 2 . Finally, the TPE correction 
changes the unpolarized e-p reduced cross section, er rc d, 
which appears as part of the denominator. This involves 
both replacing the Born form factors with the generalized 



7 



form factors, and introducing new terms related the new 
amplitude that appears when including TPE. 

We evaluate the corrections to Apy using the TPE 
calculations of Refs. [H, EH, which give identical results 
for the Q range of interest. The new terms, A' M and 
A' A , have small contributions, below 0.1% to Apy for 
Q < 1 GcV/c. While these terms become comparable in 
size to the corrections to Aa,Ae, and Am at higher Q 
values, they are negligible at low Q. The effect of the 
generalized form factors on A E , A M , and A a is larger, 
but still relatively small; on the order of 1% for Q < 
1 GeV/c. The final TPE contribution, the correction to 
cr r od (the denominator of Eq. [TJ) can be 2-3% at small e 
values. Due to cancellation between the different terms, 
the combined effect of TPE on Apy is ^1% for Q < 
1 GcV/c. 

However, the TPE correction to the Apy n is not the 
complete story; one must also take into account the in- 
direct impact of TPE corrections on the extraction of 
Ap'y . In the past, TPE contributions were neglected 
in extracting the nucleon electromagnetic form factors, 
and so the calculated value of Apy is not the correct 
value for Apy " . The TPE corrections to Ge p and Gm p 
change Apy by 1-2% for Q below 1 GeV/c, and as much 
as 5% for Q = 2 GcV/c. The largest corrections com- 
ing for small scattering angles, where the precision of 
the completed and planned measurements is the highest, 
and so one must apply TPE correction to the extraction 
of the nucleon electromagnetic form factors, as done in 
Section [TlAl 

In summary, the largest effects are the corrections to 
the Roscnbluth extracted values of Ge p and Gm p , and 
the application of the TPE effects to the denominator of 
Eq. [TJ The TPE effect on the denominator are identical 
to those in the unpolarized cross section measurements, 
and so a model-independent extraction of the denom- 
inator can be achieved by using uncorrected e-p cross 
section. We have provided TPE-corrected fits to the 
form factors, as well as uncorrected fits to the unpolar- 
ized cross sections, which allow these corrections to be 
applied without requiring an explicit calculation of the 
TPE corrections to Apy . If one neglects the remain- 
ing corrections to the numerator of Eq. [TJ the result is 
within 1% of the full calculation. However, a simple lin- 
ear parameterization of these remaining terms provides 
a calculation of Apy that is within 0.2% of the full cal- 
culation: 



A 



PV 



A 



pv 



[l + (C +ed)] 



(9) 



where C Q = .013 - .022Q, d = -.010 + .018Q, with Q 
in GcV/c. 

Figure 2] compares various approximations for Apy to 
the full calculation explicitly including TPE corrections. 
The top left panel shows that the correction to the Born 
value is small, due to the relatively small direct TPE 
contributions, and the cancellation between TPE contri- 
butions to different terms. The bottom left plot shows 
the error made when neglecting TPE corrections in both 
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FIG. 4: Comparisons of different calculations of Apy to the 
calculation including the full TPE effects. Top left plot is 
Apy~ n , bottom left is neglecting TPE in both the extraction 
of the EM form factors and in calculating Apv, i.e. the pro- 
cedure used in analyzing previous experimental results. The 
top right is the approximation presented here, neglecting the 
additional parameterization of the TPE effect on the numer- 
ator of Eq. [TJ and the bottom right is the final prescription, 
including this correction (Eq. [9j 



the calculation of Apy and the extraction of the EM 
form factors. The right hand plots show the approxi- 
mation discussed in this paper, neglecting the additional 
correction due to the effect on the numerator in Eq. [JJ 
(top figure), and including the parameterization of this 
correction from Eq. [9] (bottom figure). 



C. Extension to larger Q 

While corrections to individual terms in Eq. [JJ can be 
at the 1-2% level, and additional corrections due to TPE 
effects in the extraction of the Born EM form factors can 
be even larger, significant cancellation between different 
terms yields a total correction that is typically below 1% 
for Q < 1 GcV/c. After applying the TPE corrections 
as discussed above, the uncertainties in the TPE correc- 
tions for Q < 1 GeV/c are dominated by the uncertainty 
in extracting the TPE-corrected form factors. This un- 
certainty is taken into account in the typical 1.5% un- 
certainty assumed for radiative corrections, and thus no 
additional uncertainty need be applied. 

At larger Q, these corrections grow significantly, as 
shown in Fig. 01 and the total error made in neglecting 
TPE corrections can reach 10% by Q = 2 GeV/c. The 
procedure described here is provides a correction good 
to 0.2% up to Q = 1 GeV/c, and 1% up to 2 GeV/c. 
At higher Q, the corrections become even larger, and 
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the calculation of TPE corrections becomes less reliable. 
An estimate of the contributions from an intermediate A 
in the box diagram (5l| indicates that this contribution 
is less than 0.3% for Q 2 < 1 (GeV/c) 2 , while at Q 2 = 
3 (GeV/c) 2 , the contribution is as large as 2%, and is 
significantly more model dependent. 



IV. SUMMARY OF THE PROCEDURE 

The final prescription involves evaluating the terms in 
Eq. [U using TPE-correctcd fits for the nucleon form fac- 
tors in the terms Ae, Am, and A a, the TPE-uncorrected 
fits to CT rc d for the denominator, and applying the correc- 
tion from Eq. [5] to account for the TPE corrections for 
the terms in the numerator of Eq. [T] and the additional 
terms A' M and A' A Q. Without this final correction, 
the approximation is valid to better than 1% for Q val- 
ues from 0.3-1.0 GeV/c, and better than 0.5% except for 
Q « 1 GeV/c and e < 0.5. With this correction, the 
approximation is good to 0.2%. 

To get the overall error of the term Eq. [6l one should 
quadratically add the following contributions: 

- the effect of the error of Ge p (Fig. [T]) 

- the effect of the error of Gm p (Fig. [T]) 

- the effect of the error of G B „ (Fig. \2§ 

- the effect of the error of Gmu (1-5%) 

- the effect of the error of the e-p cross section (the 
denominator of Eqs. [Q1]) (Fig. [3]) 

- the uncertainty associated with neglected TPE correc- 
tions 

Note that in evaluating the error due to Ge p and Gm p , 
the values of the form factors are changed only in the nu- 
merator of Eq. [TJ the value of cr ro d is left unchanged, as 
it's contribution to the uncertainty is treated separately 
fSec. Ill Cj) . For the complete analysis of the uncertainly 
of PV experiments, one must of course add the uncertain- 
ties stemming from uncertainty in 0\y and G\ as well as 
uncertainty in the scattering kinematics. 

Finally, one obtains the term involving the strange 
form factors by equating the term in Eq. [5] with (Apy — 
A s py°). Thus, the uncertainty in cr re d enters again when 
isolating the linear combination Ge s + t)Gms- Because 
Apv ~ Ap\~?, t ne 1-2% overall scale uncertainty on the 
extracted value of Ge s + vGm s will always be very small 
compared to the effect of the uncertainty of a le d on Apy°, 
and so again these uncertainties can be treated as uncor- 
rected without significant effect on the final uncertain- 
ties. 



V. CONCLUSIONS 

We have evaluated the effect of TPE corrections on 
parity-violating elastic electron-proton scattering using 



the TPE exchange calculations of Refs. [T3, Thc di- 
rect effect of TPE on the parity violating asymmetry is 
small, £1% for Q < 1 GeV/c. However, the effect of TPE 
on the Rosenbluth extractions of Ge p and Gm p , which 
are needed to extract the strangeness contribution from 
the asymmetry, can be significant, and should be taken 
into account in the analysis of the parity-violating mea- 
surements. We have provided fits to thc form factors and 
their uncertainties, and provided a prescription to allows 
for an extraction of the strangeness form factors without 
explicitly requiring a calculation of the TPE exchange 
effects. As we have shown, this prescription provides an 
excellent approximation to the full procedure, based on 
tests performed using the full TPE calculation. This pro- 
vides a common set of form factors and uncertainties for 
the analysis of low-Q parity violating measurements, as 
well as a consistent application of TPE corrections on 
the extraction of the strangeness form factors. This ap- 
proach has the advantage that it can be applied without 
requiring an explicit calculation of the TPE amplitudes, 
while providing an approximation of the the TPE cor- 
rections to better than 0.2%, with model dependence in 
the TPE correction that is consistent with the assumed 
uncertainties due to RC for the measurements. 

These TPE-correctcd form factors are needed for deter- 
mining the value and uncertainty in Apy for the case of 
no strange quark contributions. They are also thc true, 
Born form factors that are related to the structure of 
the nucleon, and which should be used in thc analysis of 
other experiments. However, the effect of TPE is differ- 
ent in different observables, and one must consider if the 
Born form factors are the correct input in the case being 
considered. For example, many analyses such as Rosen- 
bluth separations in quasielastic A(e,e'p) scattering [53[ 
or the extraction of Gmu from polarization [34l | or ratio 
measurements 32[, require knowledge of the e-p cross 



section to extract information on other quantities. If one 
uses the TPE-corrected form factors, then one must in- 
clude TPE corrections in calculating the cross section. In 
such cases, it is simpler and less mo del- dependent to use 
the fits to the uncorrected cross sections. Other cases, 
such as quasielastic neutrino scattering or the case of 
parity-violating electron scattering considered here, will 
not have the same TPE effects and need to be evaluated 
with care. 
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